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The final inspiral phase in the evolution of a compact binary consisting of black holes and/or 
neutron stars is among the most probable events that a network of ground-based interferometric 
gravitational wave detectors is likely to observe. Gravitational radiation emitted during this phase 
will have to be dug out of noise by matched-filtering (correlating) the detector output with a bank 
of several 10 5 templates, making the computational resources required quite demanding, though not 
formidable. We propose an interpolation method for evaluating the correlation between template 
waveforms and the detector output and show that the method is effective in substantially reducing 
the number of templates required. Indeed, the number of templates needed could be a factor 
^ . ~ 4 smaller than required by the usual approach, when the minimal overlap between the template 

' bank and an arbitrary signal (the so-called minimal match) is 0.97. The method is amenable to 

easy implementation, and the various detector projects might benefit by adopting it to reduce the 
computational costs of inspiraling neutron star and black hole binary search. 
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I. INTRODUCTION 
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Black hole (BH) and neutron star (NS) binaries end their lives by emitting a chirp-like burst of gravitational wave 
(GW) radiation. The inspiral radiation from NS-NS sources located within a distance of 20 Mpc, and BH-BH sources 
within 200 Mpc, should be observable by a first-generation network of ground-based long-baseline interferometric GW 
detectors |IJ . Correlating the detector output (data) with a family of expected waveforms (templates) and using the 
largest correlator as a detection statistic, provides the best detection strategy for Gaussian colored stationary noise 
5h ■ (maximum likelihood, matched filtering |2|), but as we will not know the waveform parameters beforehand, a bank of 
jpP" 10-100 thousands of templates will have to be used so as not to miss out any source ||. Data analysis groups have 
J> , been, therefore, steadily looking for data analysis economization strategies. 

In setting up a bank of templates to search for a signal with uknown parameters the quantity of interest Q is the 
(frequency domain) scalar product between the signal h(-) and the template(s) #(•)■ The scalar product defined in 
■ this context is similar to the usual vector space definition, except that the measure is derived from the detector noise 
power spectral density S n (f): 



(h,g) =2 



wa^ + G.c. 



(i.i) 



1/2 

where C.C. denotes the complex conjugate. It is customary to use unit-norm templates, for which \\g\\ = (g, g) ' = 1. 
Under this assumption the scalar product is also called deflection, and denoted by d j^j. The normalized deflection 
obtained by further dividing (g,g) by \\h\ \ is the overlap or ambiguity function. It is a real number in [0, 1], and is a 
measure of how similar the template and the signal are. 

Given a binary inspiral signal and a template, the (nominal) times at which they coalesce and the phases at the 
times of coalescence are unimportant for the purpose of data analysis. We can, therefore, preliminarily maximise the 
overlap over these parameters . Such a maximised overlap is called the match Q , and is a function of the remaining 
(intrinsic) source and template parameters S. 

The minimal match T is the smallest overlap of a signal of arbitrary parameters with the template closest to it in 
the template bank Q. It is immediately related to the fraction (1 — T 3 ) of potentially observable sources which might 
be missed out 101 . 



In a recent paper |10 it has been pointed out that (under certain assumptions) the match is a quasi band-limited 
function of the difference between the intrinsic source and template parameters. As such, it can be approximated 
by a cardinal interpolating expansion, which uses only a fraction of the templates needed by the std. lattice for a 



2 



prescribed minimal match. The same economized template set can be used to build a cardinal-interpolated formula 
for the partially maximised correlator, under the same minimal match constraint. 

This was shown in jl(| for the simplest (but rather unrealistic) case of NewtonianfOPN) chirp signals. At the lowest 
Newtonian (OPN) order an inspiral waveform depends only on a single parameter . In this one-parameter template 
space, the cardinal-interpolation method JTof leads to a reduction in the number of templates by a factor w 1.4 for 
r = 0.97, and the resulting template density is close to an absolute minimum set by the theory of quasi-band-limited 
functions. As anticipated in p0| , one expects an even larger reduction in the number of templates for multi-parameter 
PN models, in view of well known properties of cardinal interpolating expansion in several dimensions (see references 
quoted in @). 

In this Rapid Communication we present results for the case of first-post-Newtonian (1PN) signals and templates, 
which support the above expectation. We show that in the two-dimensional 1PN template space the reduction in the 
number of templates is by a factor w 4 for T — 0.97. A similar reduction in the computational resources is implied, 
which we believe the interferometric detector projects such as TAMA, GEO, LIGO and VIRGO might benefit from. 
In the following we present only essential ideas and results, deferring the details to a longer version. Throughout this 
paper we use geometrized units (i.e., G = c = 1). 



II. THE POST-NEWTONIAN WAVEFORM 



We shall use the restricted PN approximation, where the amplitude of the waveform is kept to lowest (Newtonian) 
PN order, while the phase is expanded to the highest PN order available |Q. 

The evolution of the instantaneous orbital phase (p of an (adiabatically) inspiraling compact circular-orbit binary 
is obtained (in parametric form) by solving a pair of ordinary differential equations |13| ]: 

dt m dE dip v 3 dE (2 1) 

dv F dv ' dv F dv ' 

where v is the gauge independent relative velocity of the two stars, F{v) is the gravitational wave flux (luminosity), and 
E(v) the dimensionless relativistic binding energy of the system, related by F(v ) = —m(dE/dt). The instantaneous 



gravitational wave frequency few is related to the orbital phase, and hence via (2.1) to v, by 



JGW=n -77 = • (2.2) 

dt irm 

The flux and energy functions are presently known to order v 5 (2.5 PN) for non-spinning binaries H|. 

The gravitational wave (dominant, m = 2 multipole, 2nd harmonic of orbital frequency) emitted by an inspiraling 
compact binary and sensed by an interferometric antenna is described by the waveform |lq| : 

hit) = 4Cr)—v 2 (t) cos[2cp(t)l (2.3) 
r 

where r is the distance to the source, and C is a constant in the range [0, 1] dependent on the relative source/detector 
orientation with r.m.s value 2/5 (average over all orientations and wave polarisations) p6[ . 

It is straightforward to obtain an explicit frequency-domain (Fourier transform) representation of the waveform 



(2.3) using the stationary phase approximation (17J, hi 



h{f) = 2C V ^- ^ {tf) e <hM«/W4, (2.4) 
where tf is the stationary point Jl9| and an overdot denotes the derivative w.r.t. t. On substituting for the statio nary 



point tf, and consistently using the available PN expansions of the flux and energy functions, the phase in Eq. (2.4) 
takes on the simple form 



1>,{t f ) = 2nfT -<h + 2 ^ dv(vj-v 3 )^^ = 27rfT c -^c + J2^k(f)r k , (2.5) 

Jv f o 

where Vf = v(tf) = (irmf) 1 / 3 , Tc is the (nominal) time of coalescence po[ , <fic is the phase at t — Tc- The r k are 
the so-called (dimensionless) PN chirp times J2l|, |22| : 

5 i „ 2 5 /743 11 \ , tf , 5 /3058673 5429 617 2 \ ,„ „. 
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where vo = (irmfo) 1 ^ 3 , fo is a scaling frequency (to be specified below), and the \E , / C are given, in terms of the scaled 
frequency v = f /fo, by: 



(» 



A/3' 



(2.7) 



In the restricted PN approximation the amplitude in (|2.4| ) depends on frequency simply through a factor / 7 / 6 . 
III. THE FIRST POST-NEWTONIAN MATCH AND THE STANDARD TEMPLATE BANK 



In the stationary phase restricted PN approximation for h, the match (or reduced normalized deflection, in the 
jargon of |fTofl ), can be written: 



D = max 

AT C 



fs f-7/3 

At L i(277/AT c +A*) 



fi 



fs f~ 7 / 3 

df sM) 



(3.1) 



where [/;, / s ] is the antenna spectral window, AT C is the difference in the coalescence times of the signal and template, 
and, at the 1PN level of approximation, 
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1 



-At 2 



(3.2) 



where At 1 — t 1 s — Tj,, i — 0,2, are the differences in the chirp times of the source (S) and template (T), respectively, 
v = f/fo, and fo is the frequency at which S n (f) is minimum. In the sequel we shall adopt the LIGO form [p3| of 
the (one-sided) noise power spectral density, 



(3.3) 



The reduced deflection Eq. (3.1) depends on fi, f s and fo only through the integration limits, which become v x = fi/ fo 
and v s = fs/fo. 

Following Owen |?J we shall perform a coordinate transformation which, asymptotically in a neighbourhood of the 
peak Ar° = At 2 = 0, re-shapes the surface z = D into a convex circular paraboloid: 



At = Az cos - AAx 1 sin i?, 
At 2 = Aa: sin-d + AAx 1 cosd. 



(3.4) 



The above transformation depends on the choice of fo, v\ and v s . In the following, for illustrative purposes, we 
let f = 200 Hz (initial LIGO), v { = 0.2 and v s = 4, whereby d « 0.489 and A = 15.203. The resulting function 
D(Ax° , Ax 1 ) is shown in Fig. 1. 

In the standard lattice approach, the template spacing 8l is obtained by enforcing the minimal match condition, 
whereby for any admissible source, there exists at least one template such that D > T. The curves D — T, shown 
in the inset of Fig. 1, are nearly circular |^4| down to T = 0.95, which includes all values of practical interest for a 
single-step search. Thus, Sl is the side-length of the inscribed square. In the inset of Fig. 2 we display 8l as a function 

of r. 

IV. THE CARDINAL INTERPOLATION FOR THE 1-PN MATCH 

In the cardinal-interpolation approach the template spacing 8c is obtained by enforcing the minimal match condition 
on the cardinal expansion of the match pQ]. The 2D cardinal expansion of the 1PN match is: 



— oo. + oo 



D B (Ax°,Ax 1 ) 



DeiAx^Axl) sine 



8c 



(Ax° - Axt 



8c 



(A^-Axr 



(4.1) 
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where Ax^ +1 — Ax^ = Ax^ +1 — Ax* = 5c is the cardinal spacing to be determined. This expansion represents exactly 
the function: 

Db = V^a,-] [w (jf) ^[A^A^Aa?)} , (4.2) 



where T is the Fourier-transform operator, 



Ay0 Ay ix ( 1, |Ay°| <B,|Ay 1 | < fl, 



^If-f = (4-3) 
' 1 0, elsewhere, 



and B = (25c) 1 > Similar to the Newtonian case @, the 1PN match Eq. is a quasi-band-limited function in 

the L°° norm, namely, 

37, B c eK+ : sup |.D - 5 B | = exp[- 7 (S - B c )]. (4.4) 



The cardinal expansion Eq. (4.1) can be shown to approximate the function Eq. (13. 1) i n the L°° (and L 2 ) norm, using 
the minimum 2D sample density 1 / <& c ompatible with a prescribed accuracy (see |10f , and references cited therein) . 



The exponential decay of the error (4.4) is exemplified in Fig. 3 



We shall now discuss the template spacing to be used in Eq. (4.1), for a prescribed minimal match T 



The source and template parameters can be always conveniently written as follows: 

xs = XQ+a s Sc, x T = XQ + d T 6c, (4.5) 

where ctglp 6 [0, 1] and xq correspond to the lower left corner of a suitable template-space cell. 

For any given value of 5c, as the source coordinates xs move in the template-space cell, the location x1f ax where the 
cardinal-interpolated match Db attains its maximum value and the maximum value itself D B lax change. The minimal 
match condition should be enforced for the special value <3* of as where the interpolated maximum is a minimum. 

Numerical experiments show that a* = (0.5,0.5), which corresponds to the center of the cell, below some critical 
value 5 X of 5c- Above this critical value, a bifurcation occurs, i.e., two (equal) minima exist |p5| located on the 
descending diagonal of the cell, where 

a| = 0.5 + £, 4 = 0.5-£, (4.6) 

symmetrically with respect to the cell center, viz. at: 

t = ±&{6c), USc) e [0,0.5]. (4.7) 

This is illustrated in Fig. 4, where we show the density plots of D B ax as a function of as in the fundamental cell, for 
the representative cases 5c — 0.078 < 5 X and 5c — 0.099 > 5 y 



The values of D B lax on the descending diagonal of the cell, where the minima are located and Eq. (4.6) holds, are 
plotted in Fig. 5, as a function of the displacement £ off the cell center, for several values of 5c- The values of D B lax 
at the minima of these curves, occurring at £ = £*(5c)> are the minimal matches. 

It is seen, e.g., that for T = 0.97 one has 5c ~ 0.104. By comparison with Fig. 2, where 5l — 0.052 for T = 0.97, 
we deduce that 5c /5l ~ 2. Hence the reduction in the 2D template density is of the order of (5c/5l) 2 ~ 4, which is 
quite remarkable. It is also seen that, similar to the 0PN case, the 1PN template density reduction with respect to 
the plain lattice increases with T, as shown in Fig. 2. 



V. CONCLUSIONS 



The number of 1PN templates required to keep the minimal match above a given threshold T can be signifi- 
cantly (by a factor w 4 at T ks 0.97) reduced by using cardinal interpolation. The statistical properties of the 1PN 
cardinal-interpolated (partially maximised) correlator bank will be the subject of a forthcoming paper. Extension to 
2.5PN templates should be straightforward, in principle, using the (almost-flat parameter space manifold) coordinates 
introduced by Tagoshi and Tanaka p6| . 
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Figure Captions 



Fig. 1 - The function z = D(Ax° , Ax 2 ), and some of its contour levels (inset). 



Fig. 2 - The plain 1PN lattice spacing Sl (inset), and the 2D template density reduction (Sc/5s) 2 , achieved by use 
of cardinal interpolation, both plotted as functions of the match T. 



Fig. 3 - The L°° error as a function of 5 C 2 . 



Fig. 4 - Density plot of D max vs. a s for 5 C = 0.0078 (left) and S c = 0.0099 (right). 



Fig. 5 - Dg ax vs. £, for several values of 5c- 




Fig. 1 - The function z = D(Ax°, Ax 1 ), 
and some of its contour levels (inset). 




Fig. 2 - The plain 1PN lattice spacing #l (inset) and the 2D template 

2 

density reduction (S^/Si) achieved by use of cardinal interpolation, 
both plotted as functions of the minimal match T. 
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Fig. 5 - D B vs. £ for several values of £ c . 



